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ABSTRACT 

We present a study of ionized gas, PAHs, and H2 emission in the halos of 
three edge-on galaxies, NGC 891, NGC 5775 and NGC 3044, based on 10-20 fim 
Spitzer Space Telescope spectra. The [NeIII]/[NeII] ratio, an excellent measure 
of radiation hardness, rises with z in the halo of NGC 891. It is also higher in 
the halo of NGC 5775 than in the disk. NGC 3044 presents a more confusing 
situation. To explain the [Ne III]/[Ne II] as well as optical line ratio behavior in 
NGC 891, we carry out a simple exploration of parameter space with CLOUDY, 
which indicates a large increase in radiation temperature with height. Illustrative 
examples of physical models using a Monte Carlo radiative transfer code show 
that the rising neon ratio may be explained by adding a vertically extended, hot 
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stellar source to a thin disk of massive stars. However, several other sources of 
hard spectra may be relevant. PAH features have scale heights of 430-530 pc 
in NGC 891 and 720-1080 pc in NGC 5775, suggesting they can be transported 
by disk-halo flows. Within NGC 891 and NGC 5775, scale heights are similar 
for all PAHs. For NGC 891, the scale heights exceed that of 8 /xm emission, 
indicating a transition from more ionized to more neutral PAHs with height. 
Most PAH equivalent widths are higher in the halos. H2 17.03 /im emission with 
scale heights of 550 — 580 pc in NGC 891 and 850 pc in NGC 5775 suggests a 
molecular component in a surprisingly thick layer. 

Subject headings: galaxies: ISM — galaxies: spiral — galaxies: individual (NGC 
891, NGC 5775, NGC 3044) — methods: numerical 



Introduction 



Gaseous halos of spiral galaxies have grown in importance in recent years as they are 
at the interface between galaxies and their environments, and as such are the sites of many 
physical processes that dictate galaxy evolution. Supernovae and stel lar winds in disks 



can p ump mass and energy into halos setting up a disk-halo cycle of gas (INorman fc Ikeuchi 



■f 

a) 



19891 ). depending on the level of star formation activity. Halos may also con tain gas accreting 



onto galaxies, either primordial or from companions JSancisi et al.lbood l Gas from these 
two sources may interact dynamically, thermally, and chemically. Such feedback from star 
formation is believ ed to be important f or understanding the growth of galaxies and their 



current properties (jWhite fc Frenklll991 



The discovery of thick layers of interstellar gas and dust, whether referred to as "halos" 
or "extraplanar" components, has opened up a window on these issues. In the past twenty 
years or so, such thick lay ers have been discovered in just about every co mponent of the ISM, 
especially in X-rays (e.g.. Tiillmann et al.ll2006bl: [Strickland et al.ll2004[). HI fe.g.. O osterloo 



Fraternali, & Sancisi 2007), radio continuum (e.g., Dahlem. Lisenfeld. fc Rossa 2006h. diffuse 
ioniz e d gas (DIG: e.g., Irossei fc Dettmar 2003b; Rand 19961 ). and dust ( Howk fc Savage 
I999I : lAlton et al.lll998l . l2000h . Particularly for the DIG, X-ray, radio continuum and dust 



components, their brightness and extent correlate with the level of disk star formation , both 



within and among galaxies, indicating an origin in a disk-halo flow (Randl 



& Dettmar 



2003a 



Tiillmann et al. 



2006 



3 



Strickland et al.ll200 



3 



19961 



Dahlem et al. 



20o"i: Howk & 



1998a: Rossa 



bavage Il999h . Yet, some halo gas may be due to continued primordial infall onto galaxy disks, 
which some lines of ev idence indicate is necessarv to maintain star formation fe.g.. Chiappini. 



Matteucci, & Gratton 1997; Bauermeister. Blitz, fc MalbolO ). Also, the metallicities of 



many 
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Galactic High Velocity Clouds are as low as 0.1 Z© (IWakkerl 12004) , possibly indicat ing: a 



mix ing of primordial and processed gas. Finally, lagging DIG (e.g.. lHeald et al.l 120071 ) and 



HI (Oosterloo et al.l 2007 ) halos have also been interpreted as a sign ature of such ra ixing 
(IFraternali &: Binneyl l2008l ) . although other explanations are possible (IBenjami ^ l2002h . 



1.1. Spectroscopy of DIG 



The DIG is a valuable tracer of energetic processes occurring in gaseous halos, largely 
through the substantial leverage provided by observations of emission lines. Studies of halo or 
"extraplanar" DIG in external edge-ons have been r notivated by the detailed characteriz ation 
of the DIG (or Reynolds Layer) of the Milky Way jRevnoldslll99ol : iHaffner et aUboOoh . and 
the questions that have arisen therefrom. One of the most important questions is what keeps 
the DIG ionized. 



Despite the large vertical extent of extraplanar DIG layers in many galaxies (e.g.. iRand 



19961 ). there is much evidence from various diagnostic line ratios that the primary source 



of ionization is radiation leaking from a thin disk of massive stars (e.g.. iRand 1997 ^ 1998b 



Greenawalt. Walterbos. fc BraunI 119971: 



Haffner. Revnolds. fc Tuftelll999 



Collins fc Rand 



2001 



Otte. Gallagher, fc Reynolds 



2002 



Hoopes &: Walterbos 



2003 



;_Madsen, Reynolds, &; 



Haffner |2006| ). In the Milky Way, such radiation is also the only source that can meet the 
energetic requirement of keeping the DIG layer ionized ( Reynold^ 1990). Most commonly 
observed is that the ratios of [S II]AA6716, 6731 and [N II]AA6548, 6583 to Ha generally in- 
crease with distance from the suspected ionizing source, in the Milky Way and in external 
galaxy disks and halos (all emission lines, and ionization energies necessary to create the 
responsible ions, relevant for th is paper are listed in Table 1 ) . This is a sig nature of a falling 



ionization parameter, U, (e.g., iDomgorgen &: Mathislll994l : lRandlll998bl ) but other factors 



have been a rgued to affect line ratios significantly - na mely elevated gas temperatures in the 



diffuse gas ([Reynolds et al. 



20011: 



Haffner et al 



19991') and hardening of the radiation field 



during its propagation (jHoopes &: Walterbos 2003 : Wood &: Mathis 2004, hereafter WM) 



The situation is further complicated by the spatial b ehavior of other lin e ratios, such as 



[OIII]A5007/H/3, [OI]A6300/Ha, and HeIA5876/Ha fsee lHaffner et al.ll2009l for a summarvV 
The first of these ratios, in particular, shows behavior in many galaxies suggesting that some 
DIG ionization occurs through processes other than photo-ionization by a thin disk of massive 



or stays relatively constant with distance from the ionizing stars (e.g., 


Rand 


1998b. 


2000: 


Galarza. Walterbos. & Braun 


1999: lotte et al. 12002: 


Hoopes & Walterbos 


2003 


: Mil 


ler & 
n the 


Veilleux 


2003 


). An increase with height has also been found in the Reynolds Layer i 
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inner Galaxy (iMadsen &: Revnolda 12005 ). This may be partially explained by elevated gas 
temperatures, although I Collins &: Randl ( 120011 ) find this explanation insufficient for NGC 
5775 and UGC 10288. In cases such as these and the well-studied NGC 8 91, other sources 
of ion ization are indicated, with shocks receiving the most attention (e.g., I Collins fc Rand 
booih . Not all extraplanar DIG lay ers may require such a source, however. Some show 



a drop of [0III]/H/3 with height, z flTiillmann fc Dettmail I2OOOI : Miller fc Veilleuxl 120031 ). 
One of these is NGC 3044, especial l y clea r in the more westerly of the two long slit spectra 
presented by iTiillmann fc Dettmad (120001 ) . 



A further complication is the assumed abundan ces. In general, solar or "IS M" abun- 

Sembach 



danc es are assumed in photo-ionization models (e.g.. iDomgorgen fc Mathis 



1994 



et al. I2OOOI : WM), but one can expect some sensitivit y of line ratios t o abu ndances, either 
directly, or through variation in the cooling efficiency. ISembach et al.l (I2OOOI ) find that most 
line ratios vary inverselv with abundance in their photo-ionization models, although Dom- 
gorgen & Mathis (jl994j ) and WM find a more complicated relationship. Given the possibility 
that halo gas contains a mix of disk-halo cycled and primordial gas, abundance variations 
may well be important. 

Finally, one should mention the complication caused by extinction, which has been found 
to be signi ficant even at z = 1 — 2 kp c from the disk in several edge-on spirals, including 
NGC 891 (IHowk fc Savagelll999l . I2OOOI . and references therein). We crudely estimate in this 
paper that there should be many magnitudes of visual extinction in the midplanes of the 
galaxies we report on. Therefore, apart from the effect of extinction on line ratios, the depth 
along the line of sight to which we are probing optically likely varies significantly with z. 

With these various difficulties in interpreting optical line ratios, it is desirable to have a 
diagnostic of ionization with a more straight forward interpretat ion. Such an opportunity is 
provided by the Infrared Spe ctrograph (IRS; iHouck et al.l 120041 ) on board the Spitzer Space 
Telescope ( iWerner et al.ll2004l ). which allows measurement of the ratio of the 15.56 /xm [Ne III] 
(ionization potential 41.0 eV) and 12.81 yum [Nell] (ionization potential 21.6 eV) lines, as 
well as others. This ratio provides a diagnostic of the hardness of ionizing radiation that is 
relatively insensitive to extinction, gas-phase abundances, and gas temperature (being low 
excitation lines in warm gas) - three of the biggest sources of confusion for the optical lines 
[the ratio is enhanced in low-met allicity galaxies, but this is attri buted to the associated 
harder stellar radiation fields (e.g. iHunt et al.ll2010l . iHao et al.ll2009l )]. Pioneering work with 
the Infrared Space Observatory (ISO) fir st de monstrated the po wer of this diagnostic (e.g., 
Thornlev et aPboOol : Iciveon et al.lbooi and lVerma et al.ll2003[ l. 



In lRand. Wood, fc Benjamin! (120081 . hereafter RWB) we presented our first measurements 
of this ratio in the DIG halo of NGC 891 (assumed distance 9.5 Mpc), at 2; = 1 kpc on 
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opposite sides of the disk. The ratio was found to be elevated in both extraplanar pointings 
relative to the disk. 

If the DIG is ionized by radiation leaking out of the disk, then the same stars, with the 
same metallicities, are responsible for disk and halo ionization, and our result for NGC 891 
would therefore be due to propagation effects. To understand whether such a scenario is 
viable, we combined the neon ratio with optical line ratios along a long slit running through 
the IRS pointings, and mo deled line emission from the DIG l ayer using the 2D Monte-Carlo 
radiative transfer code of IWood. Mathis. fc Ercoland (l2004j ). The code is also capable of 
3D simulations. The models include the effects of radiation field hardening and explore 
various temperatures and luminosities of the radiation field, and were tailored to the diffuse 
ISM density distribution of NGC 891. As applied to optical line ratios for the Reynolds 
layer and the halo of NGC 891, WM had found that spectral hardening reduced the need 
for a non-ionizing heat source, yet it was still not possible to match all the line ratio data, 
so that additional heating and/or a secondary ionization source still seems to be required. 
With the neon ratio included, RWB found that no model could reproduce the observed low 
values of this ratio and its rise with z\ neither could they predict a rising [0III]/H/3, or the 
typical values of [O Ij/Ha, and He I/Ha. They are more successful in reproducing the rise of 
[SII]/Hq; and [Nil] /Ha, at least semi-quantitatively. We emphasize that including a non- 
ionizing heat source or a metallicity gradient should affect the optical line ratios significantly 
but not the neon ratio. These results are obviously problematic for such photo-ionization 
models. 

Clearly, one wishes to know whether the values of the neon ratio measured by RWB, and 
their enhancement off the plane, are typical of DIG layers, whether this infrared ratio reflects 
trends in optical ratios, especially [0III]/H/3 (the second ionization potential of oxygen is 
35.1 eV, nearly as high as that of neon), and whether the neon ratio in general presents such 
difficulties for photo-ionization models. We have therefore extended our observations to z = 2 
kpc on both sides of the plane in NGC 891 (where [O III]/H/3 rises more dramatically), and 
observed disk and halo positions in two other edge-ons, NGC 5775 and NGC 304 4 [we use 
the sa me distances. 24.8 Mpc and 16.1 Mpc. respectivelv. as inlCoUins et al. J2000h and Lee 



& Irwin (Il997l ))]. where coincident optical spectra along slits perpendicular to the disk ar e 



also available. The former galaxy is more actively star forming than NGC 891 (Rand 



while the long-slit spectroscopy indicates a strongly rising [OIII]/Ha with z (IRand 



19961 



2000h . 



By contrast, as mentioned above, NGC 3044 shows a falling [0III]/H/3, especially for one 
long slit whose location is the focus of our IRS observations. It also has a relatively bright 
extraplanar DIG layer (IRossa fc Dettmaiil2003al ). 



We reconsider photo-ionization models in light of our results, focusing on NGC 891. 
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We b egin with a general exploration of parameter space using the CLOUDY code f Ferland 



et al. Il998l ). asking the simple question of whether there is any combination of ionization 
parameter, radiation temperature and gas temperature that can reproduce the infrared and 
optical line ratios in NGC 891 independently at each height. That exploration opens up a 
wide range of parameter space with a new set of constraints to consider in physical models, 
and suggests a tractable way to explore multiple ionization sources in future work. Consistent 
with these findings, we present an initial exploration of this space with new Monte Carlo 
radiative transfer models that include a second, vertically extended stellar ionizing source. 



1.2. Dusty Halos 



The dust content of gaseous halos can provide cl ues to their origin. T here is plenty of 
evide nce for dusty halos th rough both emission (e.g., IWhaley et al.ll2009l ) and absorption 



Howk fc Savagelll999l ). That dusty halos are powered by star formation driven disk- 



(e.g. 

halo cycling of gas is suggested by the correla tion among galaxies b etween the presence of 



extraplanar Ha extinction and emission (e.g., iHowk fc Savagelll999l ). It seems most likely 



that dust is transported upwards w ith the gas. However, other explanations for halo dust are 
possibl e fseelHowk &: Savagelll997l for a general discussion) , including radiative acceleration of 
grains (IFranco et al.lll99ll ). which may lift clouds to heights of a few hundred pc through gas- 
dust coupling. It has also b een argued th at grain lifetimes are much shorter than the typical 
residence time in the ISM (jDrainell2003l ). implying much dust formation and modification 
in the ISM. It seems more likely that such dust formation would take place in the denser 
disk environment than in halos, although the energetic processes involved in disk-halo flows 
may modify dust properties. This latter point may be important for understanding dust 
evolution and provides motivation for studies such as ours. 

The IRS provides the ability to study emission from Polycyclic A romatic Hydrocar- 



bons fPAHs: Leger &: Puget 



Bauschlicher 2001a 



1984 



AUamandola. Tie 



Bakes et al. l2001b: Draine fc Li 



ens. &: Barkerlll985l: Bakes. Tielens, & 
20071 ) - molecules up to about 20A in 



size containing up to a few thousand C atoms which in SINGS galaxies are found to con- 
tribute a few percent of the dust mass (jPraine et al.ll2007l ). Emission features in the 10 — 20 
lim ra nge are thought to arise from more neut ral PAHs relative to features in the 6 — 9 fim 
range JPraine fc LilboO?! : Icalliano et allbood l 



It is uncertain how PAHs relate to the larger dust grains, but it has become clear in 
recent years that PAH emission is suppressed rel ative to the emission from larger grains in the 
immediate vicinity of recently formed stars fe.g.. lLebouteiller et al.ll2007), a re sult attributed 



to destruction of PAH molecules by UV radiation (e.g., iGordon et al 



2008|). A key piece 
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of evidence for this picture is the decrease of PAH Equivalent Widths (EW) with increasing 
radiation field hardn ess (at least in the 20 — 40 eV range) , as measured by [Ne III] /[Ne II], in 



star forming regions (iMadden et al.ll2006l : IWu et al.ll2006 



Engelbracht et al.ll2008[). although 



there seems to be no correlation below [Ne III]/[Ne II] 
20081 ). the regime that will be relevant here. 



1 (IBrandl et al. 



2006 



Gordon et al. 



In halos, radiation fields are much weaker, whatever their hardness, and it is less likely 
that processing of grains by radiation will be important. Furthermore, it should be pointed 
out that in the general ISM as well as in halos, radiation field intensities are low enough so 
that even the continuum in the 10 — 20 /im range should be domin ated by grain s heat ed by 
single photons rather than ones in thermal equilibrium, according to lDraine &: Lil (120071 ). We 
should therefore expect that the PAH spectrum should not vary significantly with radiation 
hardness or intensity for the kinds of environments relevant here. Rather, relative variations 
in PAH emission are likely to reflect real changes in the PAH population, and would point to 
other energetic processes likely to be relevant for halos, such as the aforementioned photo- 
levitation, or shocks, where PAHs may be produced (by grain colhsions or sputtering: Jones. 
Tielens, & Hollenbach 1996 ). as well as destroyed ( Micelotta. Jones, fc Tieleiis[2010 ). 



Imaging is beginning to reveal extended halos of PAHs in edge-on galaxies (e.g.. Whalev 



et al. 



2009 



for NGC 891. lHowk et al.ll2009l for NGC 5775). Spectroscopy reveals more detailed 
information but fewer studies exist. In RWB our limited spectra provided a first look at scale 
heights and EWs of PAH features in NGC 891. Elsewhere, the ha l o of th e starburst M82 
has also been s t udied spectroscopically, first by Engelbracht et al.l ( 120061 ). More recently, 
Galliano et al.l ( l2008l ) find that the intensity ratios of both the 6.2 and 7.7 /zm features 
relative to the 11.2 /xm feature decrease in this halo, which they interpret as a drop in the 
contribution from ionized PA Hs. Otherwise , they find the mix of PAHs in M82 is quite 
invariant. In the same galaxy, iBeirao et al.l (120081 ) find higher EWs of neutral PAHs, but 
not ionized PAHs, in the halo than in the disk, but they attribute this to a drop in the 
continuum longward of 10 fim relative to the PAH strength. 

With our larger set of IRS spectra we can start to examine these issues in the halos 
of more normal, albeit still quite actively star forming, edge-on galaxies, and do so in more 
depth than was possible in RWB. Scale heights of PAH emission can be compared with 
those of other vertically extended components to relate PAH and gaseous halos. Scale 
heights of various features and disk-halo contrasts of their EWs can be compared to search 
for modification of the PAH population with height. 
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2. Observations 

The data were taken on 2008 March 26 (NGC 5775), September 9 (NGC 891), and 2009 
January 7-8 (NGC 3044) (program ID40284; PI: R. Rand) using the staring mode of the 
IRS Short-High (SH) module on board Spitzer. A log of the observations is given in Table 2. 
The SH module is a cross-dispersed echelle spectrograph providing spectral coverage from 
9.9 to 19.6 //m with a resolving power of ~ 600. The aperture has dimensions of 4.7"xll.3". 
In staring mode the SH module nods between two pointings centered 1/3 and 2/3 of the way 
along the sht. 

Coordinates of all pointings are given in Table 2, while Figures 1, 2, and 3 show the 
locations of the apertures (including the NGC 891 apertures from RWB) on Ha images of 
the three galaxies, as well as the overlapping long shts discussed above. For NGC 891, the 
new observations were made at locations in the halo of NGC 891 centered at a height of 
approximately z = ±2 kpc (hereafter "east halo" and "west halo" pointings) on opposite 
sides of a location in the disk at 100" from the center of the galaxy. The reconstructed 
pointings indicate that the actual slit centers were at heights of 2.0 and 1.9 kpc on the cast 
and west sides, respectively. The NGC 5775 halo pointing is at 2; = 2 kpc on the SW side of 
the halo. Three locations in the disk below were observed (since a pointing NW of the disk 
position directly underlying the halo position would include the nucleus and thus perhaps 
not be representative of the disk, both additional disk pointings are toward the SE). For 
NGC 3044, our two halo pointings are at heights of 0.7 and 1.2 kpc from the disk on the NE 
and SW sides, respectively, and three disk locations were also observed. 

Because of the bright infrared background, separate sky observations are necessary for 
IRS exposures of faint emission, particularly when the emission is expected to fill the slit. 
These were performed contiguously with the target observations, using an exposure time 
equal to that of the halo observations for each galaxy. 
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Fig. 1. — Section of Ha image of NGC 891 (jRand et al.lll99Cll ). Boxes show IRS pointings, 
including in red the three pointings from RWB closest to the plane (for each pointing two 
overlapping boxes showing the two nods are drawn). The solid line shows the orientation of 
the sht for the optical emission line data discussed. 




Fig. 2. — IRS pointings for NGC 5775 overlaid on the Ha image from ICoUins et al.l (l2000l ). 
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Fig. 3. — IRS pointings for NGC 3044 overlaid on the Ha image from ICoUins et al.l (l2000l ). 
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The data were processed through versions S17.2 (NGC 5775), S18.1 (NGC 891) and 
SIS. 5 (NGC 3044) of the IRS reduction pipehne, providing the Basic Cahbrated Data (BCD) 
products. Post-BCD processing was carried out according to the Infrared Spectrograph Data 
Handbook Version 2.0. Charge accumulation can be present in long integrations, evidenced 
by steady increases in the signal, but no significant effect was found in our data. "Rogue" 
pixels were cleaned using the contributed IRSCLEAN_MASK software. For each of the AORs 
listed in Table 2 and for each nod, the 2-D spectra were averaged using sigma-clipping 
at the 3(7 level. For each galaxy, the resulting sky spectra for both nods were averaged 
and subtracted from all target observations. All remaining extreme pixels not removed by 
IRSCLEAN_MASK were replaced with averages of adjacent pixels in the sky-subtracted 
spectra. 

One-dimensional spectra were extracted (and calibrated based on standard star obser- 
vations) with the package SPICE (version 2.2), using the full aperture width. Note that, due 
to slit loss correction uncertainties, the flux calibration is only fully accurate if the source is 
uniformly bright and covers an area much larger than the slit. In practice, our sources usually 
show some gradient across the slit (at worst a factor of 3 for the NGC 3044 halo pointings; 
the halo pointings of NGC 891 also show a significant gradient at about the 30-50% level). 
However, for the most part we can argue that the flux corrections are cither minimal or 
accounted for in the scatter of our derived values. The two halo pointings for NGC 891 
yield similar fluxes, and our main conclusions hold regardless of the scatter between these 
pointings (the same is true of the data in RWB). For NGC 5775, the three disk pointings 
also show similar fluxes, especially for the PAH features, which have a typical dispersion 
of 5% from pointing to pointing. This value may serve as a reasonable upper limit for slit 
loss correction variations given the contribution from real structure in the disk and other 
sources of error. We therefore assume fluxes from our halo pointing are also not significantly 
affected by slit loss uncertainties. For NGC 3044, apart from such uncertainties, our results 
are unfortunately compromised by what we assume is a large degree of real variation in the 
disk pointings. 

A final spectrum for each pointing was formed by combining the resulting 1-D spectra 
for each AOR and nod. As the spectral orders overlap, the ends of each order were removed 
from the final spectra. The spectra are shown in Figures 4, 5, and 6. 
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Fig. 4. — IRS SH spectra of the east and west halo pointings of NGC 891 at 2; = 2 kpc. 
Detected gas and the five strongest PAH features are indicated in the top paneL 
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Fig. 5. — Same as Figure 4 but for the disk and halo pointings of NGC 5775. 
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Same as Figure 4 but for the disk and halo pointings of NGC 3044. 
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As in RWB, we checked that the fluxes in our halo pointings cannot be due to disk 
emission convolved with the PSF of the IRS SH module using the Sj^zteer contributed software 
STINYTIM. In all cases, such a contribution from disk emission to the halo fluxes is negligible 
compared to the measured values. 

Line intensities and EWs were measured with the IRAF^ program splot by summing 
pixels over the line extent and subtracting a linear baseline. We compute error bars on inten- 
sities and EWs following RWB, taking into account noise in the spectra and the estimated 
uncertainties in the flux scale of IRS data, as described in the Infrared Spectrograph Data 
Handbook Version 3.1. Random noises measured around the relatively flat continuum at 15 
Hm wavelength range from about 1 mJy for the disk pointings to 0.07 mJy for the NGC 891 
halo pointings. 

We identify the following stronger PAH features and give in parentheses the wavelength 
range used to estimate their fluxes and and to deflne a linear continuum for the EWs: 11.2 
/im (11.0-11.8 /im), 12.0 /im (11.9-12.2 /xm), 12.7 /im (12.4-13.2 /im), 16.5 /im (16.3-16.8 
yum), and 17.4 /im (17.3 — 17.7 /im). Weaker features are also seen in som e cases at 13.6, 14. 3 
and 15.9 /im. All of these features are commonly seen in galaxies (e.g., ISmith et al.l 120071 ). 
The spectra from RWB were also reanalyzed to ensure that all features were measured in a 
consistent way. 

Defining the continuum for these PAH features is difficult and in some cases the "con- 
tinuum" values may include blended emission from discrete features. It may therefore be 
instructive to consider an EW which minimizes this possible contamination. Hence, we also 
consider EWs of the combined 11.2, 12.0 and 12.7 /im and weaker 13.6 and 14.3 /xm com- 
plexes using a linear continuum from 10.8 to 14.5 /xm, and the combined 17 /xm complex 
using a linear continuum from 16.3 to 17.7 /im. Spectral decompositions of low-resolution 
IRS spectra suggest that the emission in t hese two wavelength ranges is dominated by dis- 
crete PAH features (e.g.. lSmith et al.ll2007l ). while the ranges used for continuum are indeed 



essentially free of discrete features. Even for this type of EW, however, there is concern that 
the 10.8 um continuum mav be lowered bv silicate absorption centered at 9.7 um fe.g.. Chiar 



& Tielens l2006l ). We will discuss the effects of extinction in general on our results below. 
Just as impor tantly, even the "con tinuum" in the SH band is probably due to single-photon 
PAH heating: iDraine et al.l (120071 ) conclude this to be the case for 16 /im emission observed 
in SINGS galaxies with the "blue peak-up" array of the IRS. 



^IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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3. Results 

3.1. Ionized Gas Phase Emission Lines 

We detect [Nell] 12.81 //m, [Nelll] 15.56 //m, and [SIII] 18.71 //m in all pointings. Line 
intensities and ratios are presented in Tables 3 — 9, along with upper limits on [SIV] 10.51 
/im emission. 

The data for NGC 891 also warrant a graphical presentation. [Nelll], [Nell] and [SIII] 
intensities for NGC 891 from the current data and RWB are plotted in Figure 7a while Figure 
7b shows [Ne III] /[Nell] and [S III] /[Nell] ratios. The falloff of intensity with z is not weU 
described by an exponential for some lines. Nevertheless, we report equivalent exponential 
emission scale heights in Table 10 and plot them in Figure 8 (error bars reflect the scatter in 
the data on the two sides of the disk) to give a sense of the vertical extent of the emission. 
The key result is that the trend of rising [Ne III]/[Ne II] for this galaxy found by RWB is 
continued by the z = 2 kpc pointings. The near constancy of [S III]/[Ne II] reported by 
RWB, however, is only continued in the east halo pointing at z — 2 kpc; the value in the 
west halo is significantly higher. This is somewhat surprising given the similar ionization 
potentials of S++ and Ne+. 

For NGC 5775, there is a small but significant (at about the 5a level) [Ne III] /[Nell] 
contrast between the three (remarkably consistent) disk and one halo pointings. The equiv- 
alent exponential emission scale heights for the three lines are listed in Table 10 and plotted 
in Figure 8 (error bars refiect the scatter in the disk values). For NGC 3044, the disk and 
halo pointings show a tremendous amount of scatter for the intensities and both line ratios, 
precluding any trends from being seen. 
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Fig. 7. — (a) The log of intensity (in units of 10"^'' erg cm~^ s~^ arcsec"^) vs. z for the [Ne III] 
(red open circles), [Nell] (black filled circles) and [S III] (green crosses) lines in NGC 891. 
Vertical error bars are smaller than the symbol sizes and are omitted for clarity but listed 
in Table 3. (6) [Ne III]/[Ne II] (filled black circles) and [S III]/[Ne II] (open red triangles) vs. 
z in NGC 891. Horizontal error bars on the [Nell] points in (a) indicate the slit extent in 
the direction perpendicular to the disk. 
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Fig. 8. — (a) Scale heights of gas emission hnes and PAH features in the east (circles) and 
west (squares) halos of NGC 891 and the halo of NGC 5775 (triangles), ordered by the 
wavelength of the feature. Solid symbols represent values uncorrected for extinction and 
include error bars, while open symbols include a rough extinction correction (see text). The 
NGC 891 points have been offset shghtly in wavelength for clarity. 
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As in RWB, we crudely estimate the effects of extinction at the positions of our midplane 
pointings on the above results based on CO and HI observations and assuming Galactic ex- 
tinction properties. There are many assumptions in this calculation, including the conversion 
from CO intensity to H2 column density, the gas-to-dust ratio, and the near-IR extinction 
law. Resolution mismatches between these maps and the IRS SH observations are also a 
source of error. For NGC 891, we use the same gas column densities as i n RWB. For NGC 
5775, we use SEST CO 1-0 and uniformly- weighted VLA HI maps from iLee et al.l ( 1200 ll ). 
with the caveat that the CO map has a relatively low angular resolution of 43" (the resolu- 
tion of the HI map is 13 .5"). For NGC 3044, we use the uniformly weighted VLA HI map 
from iLee &: IrwinI (119971') a t abo ut 21" resolution, and (kindly provided to us by J. Irwin) 
the CO 2-1 map from iLed (119981 ) at 21" resolution, along with an average CO 2-1/1-0 ratio 
of 0.41. We convert CO intensities to H2 colur nn densities assumin g a conversion factor 
of 3 X 10^0 mol cm-2 (K km s-^)-\ as used bv IScoville et all Jl993h for NGC 891. With 
these assumptions, we estimate total gas column densities at the position of our central disk 
pointing in each galaxy. 

From these column densities, we estimate the midplane extinction at a representative 
wavelength of 15 /zm, makin g the same assumptions as in R WB: the Galactic relation between 
gas column density and Ay ( Bohlin. Savage, fc Drakeiri978l ). the standard Galactic extinction 
law to relate Ay to Ak (e.g.. 



Draine 



I2OO3D. and an a pprox imate relatio i i A-jr^/ Afc ~ 0.4 from 
Milky W ay observations bv iRosenthal et al.l (I2OOOI ) and Ijiang et al.l (120061 ) [although the 
model of IWeingartner fc Draind (I2OOII ) for Ry = 3.1 indicates a value closer to 0.2]. We 



also assume that the dust and emitting gas are well mixed and thus only half of the dust 
lies in front of the gas on average. For NGC 891, we estimate a midplane extinction at 
15 /im of 0.5 mag, which supercedes the value reported by RWB. For NGC 5775 we find a 
value of only 0.12 mag, and for NGC 3044 the extinction is negligible. Since these values 
are subject to suc h large uncertainties, t hey can be little more than illustrative. For the 
extinction curve of [Rosenthal et al.l (l2000l ). extinction should be about 50% higher at 17 /xm, 
and perhaps 100% higher at 11 — 12 fim, due to silicate absorption at 9.7 and 18 /xm. These 
midplane extinctions are also likely to be lower limits given the low resolution of many of 
the CO and HI observations. It may well be that the bulk of the molecular and atomic gas 
is in unresolved layers for many of these maps. However, given the lack of observations of 
sufficiently high resolution, we do not attempt to quantify this possibility here by making 
any assumptions. 

A midplane extinction of 0.5 mag, more appropriate for NGC 891, reduces scale heights 
by 30%, while a value of 0.1 mag, more appropriate for NGC 5775, reduces scale heights by 
8%. For the aforementioned infrared extinction curves, extinction for the [Nell] and [Nelll] 
lines should be comparable and therefore line ratios should not be signficantly affected. 
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Extinction for the [S III] line may be somewhat higher. 

3.2. PAH Emission Features 

Tables 3, 5, 7 and 8 show the intensities and EWs of PAH features for the three galaxies, 
as well as the intensity integrated over the most feature-free spectral range, 14.5 — 15.5 /im, 
which we refer to as the 15 /im continuum. Intensities for the 11.2 /im feature include the 
secondary peak at about 11.0 //m. EWs of the combined features discussed in §2 are also 

listed, and are hereafter referred to as EWshort and EWiong- For the NGC 891 pointings in 
RWB, EW.hort = 3.64±0.11, 3.75±0.11 and 3.83±0.11, EWiong = 0.56±0.02, 0.78±0.02 and 
0.73 ±0.02, and the 15 fim continuum intensities are 800 ±61, 32 ±2, and 35 ±3 (all in units 
of 10^^^ erg cm^^ s^^ arcsec^) for the disk, cast halo at 2; = 1 kpc and west halo at z = 1 
kpc, respectively. The spectra allow us to characterize the vertical distribution of PAH and 
continuum emission and to search for general differences of PAH emission properties between 
disks and halos. 

The rather large scatter in PAH feature brightness in the disk pointings of NGC 3044 
preclude meaningful statements about the vertical distribution of emission. Hence we focus 
on NGC 891 (where the GTO spectra discussed by RWB confirm that our disk pointing is 
representative) and NGC 5775, where the three disk pointings show a very small scatter in 
PAH emission. Figure 9 shows PAH and 15 /im continuum intensities vs. 2; on a logarithmic 
scale for NGC 891. Most features can reasonably be described by exponential distributions; 
and results of exponential fits to the east and west sides separately are presented in Table 
10 and Figure 8. Scale heights range from 430 to 530 pc for PAH features (with an average 
of 475 ± 30 pc, the error accounting for the dispersion of the ten values in Table 10), similar 
to the range reported in RWB, while the 15 fim continuum has scale heights of 440 pc and 
380 pc for the east and west sides, respectively. Hence, we have a hint that PAH emission is 
more vertically extended than the continuum emission, although the uncertainties are large. 
The EW analysis below will strengthen this conclusion. Table 10 and Figure 8 also show 
scale heights for NGC 5775, assuming an exponential distribution. For these, intensities for 
the three disk positions have been averaged, and the error bars on the scale heights reflect 
the scatter and the average uncertainty of the three measurements. Here, the scale heights 
are much larger than in NGC 891, ranging from 720 to 1080 pc. But within each galaxy, 
PAH scale heights are similar, given the error bars. 
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Fig. 9. — The log of intensity (in units of 10^^^ erg cm~^ s"^ arcsec"^) vs. z for the 11.2 
(filled black circles), 12.0 (open red circles), 12.7 (blue crosses), 16.5 (cyan squares) and 17.4 
/im (magenta triangles) PAH features and the 15 /zm continuum (green diamonds) in NGC 
891. Horizontal error bars on the 11.2 /xm PAH points indicate the slit width. 
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Extinction will lower scale heights as described above for the gas lines. Scale heights of 
the 11.2 and 17.4 fim features are likely to be reduced more by extinction than those of the 
other features. 

The PAH emission scale heights in NGC 891 with no extincti on correction a r e larg er than 



that of the CO-emittin g gas layer, which has a FWHM of 225 pc fIScoville et all (Il993[ )l The 



HI layer is modeled by lOosterloo. Fraternali. fc Sancisil (120071 ) as thin and thick exponential 



layers with scale heights increasing with radius from < 0.3 to 0.5 kpc for the former and 
1.25 to 2.5 kpc for the latter, which contains 30% of the gas. Hence, the PAH emission scale 
heights are clearly larger than the CO scale height and are perhaps more comp arable t o that 



of the thin HI layer, with or without our extinction correction. For NGC 5775, llrwini (Il994j ) 
finds an exponential scale height of 1.1 kpc for the HI, while the molecular layer thickness is 
presumably substantially less but is not well known. So again, it is probably the case that 
the PAH scale heights are comparable to or somewhat less than that of the HI layer. 

More significant results can be extracted from the EWs. These are not sensitive to 
any possible extinction contrast between disk and halo, although individual values may be 
affected by the shape of the extinction curve. Also, for NGC 3044, while intensities show 
much scatter among the disk positions, the EWs are more consistent. Figure 10 shows EWs 
of the five bright PAH features, and EWshort and EWiong, plotted as a function of z for 
NGC 891, NGC 5775 and NGC 3044, respectively. These figures make it clear that, for 
most features (all five in Figure 10a, three of five in Figures 10c and lOe, and all in Figures 
10b, lOd and lOf), EWs are higher in the halo than in the disk, no matter which way they 
are measured. Note that for NGC 891 this conclusion differs from that in RWB, where a 
significant disk-halo EW contrast was found only for the 17.4 /im feature. We attribute this 
difference to the more careful EW measurements in this work. Given the scatter among 
galaxies, no feature stands out here as having a significantly higher disk-halo EW contrast 
than others. 
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Fig. 10. — Equivalent Width vs. z for PAH features, (a) 11.2 (filled black circles), 12.0 
(open red circles), 12.7 (blue crosses), 16.5 (cyan squares) and 17.4 /im (magenta triangles) 
features in NGC 891. (b) EWshort (filled black circles) and EWiong (blue crosses) in NGC 
891. (c) and (d) same as (a) and (6), respectively but for NGC 5775. (e) and (/) same as (a) 
and (b), respectively, but for NGC 3044. Horizontal error bars on the 11.2 fim PAH points 
indicate the slit extent in the direction perpendicular to the disk. 
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3.3. H2 Emission 



The H2 J = 3 — 1 17.03 fim line is detected in all pointings. H2 mid-IR emission 
lines arise f rom rotational exc itation in warm gas. Such emissi on is usually attrib uted to 
PDRs (e.g.. Ivan Dishoeckll2004j ) but may also arise in shocks (e.g.. iNeufeld et al.ll2006l ). Scale 
heights for best-fit exponentials are given in Table 10 and plotted in Figure 8. For NGC 
891, the values of 550 and 580 pc are somewhat higher than found by RWB without the 
z = 2 kpc data points. Again, the scale height in NGC 5775, 850 pc, is significantly larger. 
Extinction (see above) may lower these values somewhat. Like the PAH features, these scale 
heights are comparable to the HI scale heights in both galaxies. It is difficult to interpret 
these results in terms of column densities of infrared-emitting molecular hydrogen from only 
one spectral line. A full analysis would require several rotational transitions to be observed 



so that excitation temperatures and ortho-par a ratios could be modeled (e.g., iBrunner et al. 



20081 : iRoussel et al.l 120071 : iNeufeld et al.ll2006l ). Nevertheless, the discovery of emission from 
H2 with scale heights in the range 500 — 800 pc suggests a molecular gas component with a 
surprisingly large vertical extent. We do note, however, that the re have been indications of 



a faint, vertically ext ended CO emitting component in NGC 891 ( iGarcia-Burillo et al.lll992 



Sofue fc Nakailll993l ). Again the greater scale height in NGC 5775 suggests that disk-halo 
flows are somehow responsible for this vertical extent. 



4. DIG Photo-ionization modeling 

As stated in §1, the [Ne III]/[Ne II] ratio provides an excellent diagnostic of ionization 
as it is relatively free of effects of extinction, gas temperature and gas abundances that can 
complicate interpretation of common optical line ratios. The extraplanar pointings at z = 2 
kpc on either side of the plane of NGC 891 continue the trend of rising [Ne III]/[Ne II] with 
z found by RWB. There is a significant but much smaller rise in NGC 5775, while scatter 
in the ratio prevents any trend from being observed in NGC 3044. All values o f the r atio in 



DIG range from about 0.1 to 1. A similar trend is seen in M82 by lBeirao et al.l (l2008l ). with 
values increasing from about 0.15 in the disk to a high of 0.27 a few arcseconds above the 
disk. 

These results indicate a higher ionization state in the halos of t hese two galaxies, con - 



firming what had been suggested by the more problematic [O HI] line ( ICoUins fc Randll200ll ). 
For what follows, we focus on NGC 891, where we have by far the most information on halo 
line ratios. In addition to [Ne III]/[Ne II] values. Figure 11 shows the dependence with height 
of [NII]A6583/Ha, [S II]A6716/Ha, [O I]A6300/Ha, [O III]A5007/H/3, and He I A5876/Ha (a 
relatively clean optical diagnostic of ionization, for which extinction is the major concern 
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and the faintness of the hehum hne the major challenge) along a slit runni ng through n early 
all of our IRS pointings (as shown in Figure 1). The optical data are from iRandl (119971 ) and 



Randl fll998bf ). 
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Fig. 11.— Dependence of (a) [Ne III] /[Nell], (b) [N II]A6583/Ha, (c) [S II]A6716/Ha, (d) 
[OI]A6300/Ha, (e) [O III]A5007/H/3, and (f) He I A5876/Ha on z for the NGC 891 lo ng slit 
spectru m discussed in the text. The Ha profile, normalized to unit intensity, from iRand 
( ll998bl ). is shown in (b). Horizontal error bars in (a) and (f) reflect the extent over which 
the data have been averaged. In the other panels, the intensities were averaged over 317 pc. 
The data are from iRandl (119971 ) and iRandl ( 1l998bl ). 
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We wish to understand what sources of ionization can explain these behaviors, with 
particular attention paid to the constraints provided by the neon ratio. 

Our main focus in this paper is on stellar photo-ionization, but by necessity we attempt 
to broaden the parameter space heretofore considered in ioniz ation models in gene ral. We 



take two approaches. We begin by using the CLOUDY code (IFerland et al.l 119981 ) to ask, 
independently at each height of 0, 1 and 2 kpc, whether there are any combinations of 
radiation temperature, ionization parameter, and gas temperature that can reproduce the 
observed ratios (§4.1). By simply exploring this parameter space, we temporarily set aside 
the issue of whether any solution corresponds to a reasonable physical or "structural" model 
- that is, we do not try to explain why there might be changes in radiation temperature with 
z, or what balance of heating and cooling may lead to a certain gas temperature. However, 
the difficulty encountered in previous modeling, even before the strong constraints provided 
by the infrared lines were available, suggests that this naive approach is a useful way to 
constrain what kinds of physical models may be viable. 

Before turning to physical models, we make several remarks about shocks (§4.2). Al- 
though standard shock models in combination with photo-ionization have been explored to 
explain the optical line ratio behavior in the DIG of several edge-ons (see §1), we point out 
the difficulty in applying such models to halo environments, and that our approach offers a 
more feasible way of considering them. 

In considering secondary sources of ionization, we have been partly motivated by the 
possibility of a second, hot, vertically extended source of photo-ionization as a way of pro- 
ducing a rising [Ne III]/[Ne II] ratio with z. Indeed, our CLOUDY results will give us further 
motivation to consider such a source. As an illustration of such a physical model, in our 
se cond approach (§4.3) w e will then turn to the Monte Carlo photo-ionization code described 



m 



Wood fc Mathid (120041 ) and employed in RWB, adding such a component to the thin layer 
of massive stars we considered in that paper. This code is well suited to exploration of 
additional stellar sources. Our CLOUDY results will open up a broad range of parameter 
space to consider in such models, and it is beyond the scope of this paper to explore fully 
this parameter space - hence we restrict ourselves to a few illustrative cases. 



4.1. Line Ratio Constraints on Physical Conditions vs. Height in NGC 891 

The increasing level of ionization in the gaseous halos indicates an increase in the ion- 
ization parameter, the hardness of the radiation field, the gas temperature, or a combination 
of all three. One difficulty in explaining the observed emission lines lies in the fact that at 
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any given height there is a range of densities. There is also quite probably a range of tem- 
peratures, either due to variation in photoionization heating or the presence of other forms 
of dissipational heating, such as shocks, turbulence, cosmic rays, or magnetic reconncction, 
that are expected to become important at low densities. It is therefore to be expected that 
it will be difficult to match the observed ratios with a single model. Uniform density mod- 
els seem likely to fail, and models with a priori density variations would have to guess the 
right combinations of density variations to match all the line ratios. Past attempts at pho- 
toionization modeling have confirmed that it is challenging to match the observed emission 
fines. 

With several measured line ratios, particularly in the halo of NGC 891, this problem 
is overdetermined. In this subsection, we use these line ratios to try diagnose the local 
physical conditions in the gas as a function of height, independent of any model to produce 
these physical conditions. In an ideal situation, in which one density /temperature phase 
dominates the observed emission, one set of physical parameters would match all the line 
ratios simultaneously. Although some line ratios remain problematic, we find that four key 
line ratios can be matched simultaneously and yield a change of physical conditions versus 
height that is physically plausible. 

We have used the photoionization code CLOUDY v8.0 for a solar metallicity gas to 
calculate how measured line ratios in the gaseous halo of NGC 891 vary as a function of 
three parameters: the ionization parameter, U = ^H/icnn), the "radiation temperature" 
(for convenience we parametrize the spectral shape with a blackbody spectrum), Tf, , and the 
gas temperature, T. The ionization parameter is defined in terms of the plane-parallel fiux of 
photons with energy greater than 13.6 eV, the speed of light, c, and the particle density 
of hydrogen atoms, uh- The ionization parameter and blackbody temperature determine 
the photoionization rate of the different ions by specifying the flux and spectral shape of the 
radiation fleld. An absorbed stellar flux or spectra from X-ray emitting gas could also be 
used. 

Figure 12 shows how the [Ne III] /[Nell] ratio varies as a function of these three param- 
eters. In the range of parameters space that we considered, this ratio spans flve orders of 
magnitude! With the relatively small uncertainties in the measured line ratio at each of the 
three heights, a thin sheet of the three-dimensional parameter space is allowed. To go from 
the ratio observed at 2; = kpc to that observed at z = 2 kpc, holding two of the three 
parameters fixed, would require a whopping 1.25 dcx (xl8) increase in the gas tempera- 
ture, a 0.7 dex (x5) increase in ionization parameter, or a modest 0.1 dex (xl.25) increase 
in blackbody temperature. The relative insensitivity of this line ratio to gas temperature, 
which only occurs because of the temperature dependence of the recombination rates, makes 
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Fig. 12. — Predicted emission line ratios for [Ne III] /[Nell] as a function of gas tempera- 
ture, logT, and ionization parameter, logf/, for four values of the blackbody temperature 
logTfe =4.55, 4.65, 4.70, and 4.75. (Note that the Tb "shces" shown in each of four quadrants 
will differ in subsequent plots.) The measured ratios at 2; = 0, 1, and 2 kpc are indicated as 
diagonal dashed lines, with the solid bands indicating the range allowed when uncertainties 
are included. The z = kpc and z = 2 kpc sheets are separated by 1.25 dex in logT, 0.7 
dex in logU, and 0.1 dex in logT^. The allowed parameter space range after incorporating 
selected other line ratios is shown (by black boxes) for convenience and discussed further in 
the text and subsequent figures. 
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it a key diagnostic. 

Incorporation of the other observed gas hne ratios produces equivalent shces through 
parameter space, most of which have a stronger gas temperature dependence. The allowed 
parameter space based on measured line ratios at 2; = 0, 1, and 2 kpc are shown in Figures IS- 
IS. The ratios [SII]/Ha and [NII]/Ha behave similarly, requiring a higher gas temperature 
to match the observed data as the ionization parameter increases. [0III]/H/3, on the other 
hand, has the opposite behavior, requiring lower gas temperature as the ionization parameter 
increases. The intersection point of these three ratios already identifies a very limited section 
of parameter space. With the addition of the infrared [NeIII]/[NeII] ratio, the problem 
becomes over-determined. The fact that the parameter space allowed by the intersection 
point of these optical line ratios (as well as He I/Ha at z = 1 kpc, although the uncertainties 
are relatively large) also agrees with the [Ne III]/[Ne II] is suggestive that one type of gas, 
with a unique set of physical conditions, dominates the observed emission. (We will address 
the remaining ratios shortly.) 

Although the intersection of these four line ratios in parameter space is not perfect for 
z = kpc, it is satisfactory given the uncertainties in the abundances for NGC 891. And 
for z = 1 and 2 kpc, it is possible to identify a unique set of parameters that simultaneously 
match all four line ratios (and He I/Ha aX z = 1 kpc). The derived parameters, as a function 
of height, are given in Table 11 (the la uncertainties define a volume in parameter space - 
the quoted uncertainties are those projected onto each axis). This table also gives the derived 
value for the (plane-parallel) flux of hydrogen ionizing photons, = cuhU, normalized to a 
density of 10~^ cm~^. Since these values are derived solely from the line ratios, it is possible 
that non-physical values could have been obtained, and yet they seem plausible. First they 
indicate a drop in ionization parameter by a factor of five over a distance of 2 kpc (despite 
the increasing level of ionization). When converted to a flux of hydrogen ionizing photons, 
it indicates an incident flux at 2; = kpc of = Ti-i(2.40 ± 0.55) x 10^ photons s~^ cm~^, 
where n_i = n/^/O.lcm"^. If we assume a constant density, only 41% and 20% of the 
incident flux reaches z = 1 kpc and 2 kpc, respectively. Using a representative gas density 
scale height of = 2.S kpc (see RWB) and a fixed filling factor, these decrease to 27% and 
8% of the incident flux at z = 1 kpc and 2 kpc, respectively. Both the total photon flux 
and the order of magnitude attenuation are comparable to what is observed in the Galactic 
warm ionized medium, which has a total flux of hy drogen-ionizing photons of = 2.6 x 



10^ photons s ^ cm ^ at z = kpc (lRevnoldslll990l ) and $h ~ 0.2 x 10^ photons s ^ cm ^ 



well above the pla ne, based on the measurement of Ha emission from high velocity clouds 



( iTufte et al.ll2002l ). The values obtained for the Galaxy were obtained by measuring the Ha 
flux and not from any line ratio information. It is encouraging that two different techniques 
yield such similar numbers for two galaxies. 
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Fig. 13. — Allowed parameter space (logT and logf/) for four values of the blackbody tem- 
perature logTft =4.50, 4.55, 4.60, and 4.65. The line ratios used to constrain the parameter 
space are indicated by the dotted lines, with uncertainty ranges shown with the solid lines. 
Each line ratio is indicated by a different color with the observed value (and uncertainty) 
given to the right of the plot (for doublets, the same line is plotted as in Figure 11). A 
successful model would have all the bands (or in the three dimensional parameter space, 
sheets) intersect in a single locus. In this and subsequent plots, the line ratios [NeIII]/[Ne 
II], [S II]/Ha, [N II]/Ha, and [O III]/H/3 have a common intersection point noted with a black 
box. Uncertainties of la on the line ratios have been used to define the uncertainties on the 
parameters. The parameters for this intersection point are given in Table 11. 
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Fig. 14. — Same as Figure 13 but for the observed NGC 891 line ratios at ^ = 1 kpc, and 
for logTfe =4.60, 4.65, 4.70, and 4.75. 
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Fig. 15. — Same as Figure 14, but for observed NGC 891 line ratios at z = 2 kpc. 
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Over the same distance, the derived gas temperature in the halo of NGC 891 increases 
by 20%, from 7130 ±820 K to 8600 ±300 K. This modest increase is also consistent with the 
increas e in gas temperature that has been claimed for the Milkv Wav bv Revnolds. Haffner, 
& Tufte ( 119991 ) . who pointed out that this might be evidence for a non-photoionization source 
of heating at low densities. One test of this temperature gradient in edge-on systems would 
be a detection of the faint [Nil] 5755 A line , as the ratio of this to the brighter [Nil] line 
is an unambiguous temperature diagnostic (IReynolds et al.ll200ll ). So far, only 2a upper 
limits on the temperature in the DIG halo of NGC 891 of 13,000 K and 10,000 K on the eas t 
and west sides, respectively, have been set from the very sensitive spectrum of iRandl ( ll997l ). 
However, another wa y of estimating the gas tem perature, based on the [N II]A6583/Hq; ratio 
( IHaffner et al.lll999l ). and applied to NGC 891 ( ICoUins k. Randll200ll ). also suggests a rise 
in gas temperature consistent with the derived rise here within the errors. 

The most surprising result of our parametrized approach to constraining the physical 
conditions is the required hardening of the radiation field. Parametrized as a blackbody 
temperature, this is a 58% increase, from 36,000 K to 57,000 K. This change cannot be 
explained by the radiation hardening expected as the photoionizing flux from the base of the 
halo loses lower energy photons as the radiation percolates upward (RWB). Some additional 
source of ionizing radiation with a hard spectrum is needed. However, this source of radiation 
must not produce too much thermal heating of the gas, nor may it exceed the constraints 
on the ionization parameter. In §4.3, we explore one such possible source. 

It is fortunate that this approach, with three parameters, can match four or five line 
ratios simultaneously, but what about the ratios we have neglected so far? Of the three 
remaining ratios, the most challenging problem is to match the rising [01] /Ha at the same 
time that [O III]/H/3 is also rising. In general, [O I] emission is thought to trace of the neutral 
phase. The neutral oxygen fraction is "locked" to the neutral hydrogen fraction by charge 
exchange, but the gas must also be hot enough to excite the optical line emission. For 
that reason, it has generally been assumed that [O I] should tra ce regions of low or p artial 
ionization, in particular the warm boundaries of neutral clouds ([Reynolds et al.lll998l ). 



Figures 13-15 make it clear that to match both this line ratio and [SHI] /[Nell], while still 
maintaining the same values for [SII]/Hq; and [NII]/Ha, requires an ionization parameter 
of \ogU = —4.6 {z = kpc) or logf/ = —4.8 {z = 1 kpc). If we were to assume a two-phase 
medium with identical gas temperatures and irradiating spectral shape, this would require a 
density ratio of 32 and 20 at 2; = kpc and 1 kpc, respectively. 



contrasts for a two-phase interstellar medium (IWolfire et al. 



lese are plausible density 
20031). Howe ver, this assumes 



that the Ha and [O I] emission regions are co-spatial. [Reynolds et al.l ( 119981 ) argue that in the 
Galactic warm ionized medium most of the Ha emission arises in a different physical region 
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than the [01] emission. This would increase the [O I] /Ha ratio in the [01] emitting zone, 
requiring a different set of parameters. (For a fixed temperature and radiation temperature, 
this would yield an even lower ionization parameter.) 

In a future study we will see if a two-phase or multi-phase combination can be con- 
structed to match all of the emission line ratios. This will include the line ratios considered 
here as well as the increasing [Oil] 3727 A/H/3 emission observed in edge-on spirals f Otte 
et al. I2OOII ) - includin g along the same slit loc ation in NGC 891 as considered here - as well 



as the Galactic WIM (IMierkiewicz et all 120061 ). 



4.2. On Shocks 



It has also b een suggested that shocks may play a role in increasing the i onization level 
with height (e.g. lRandlll998bl IColhns fc RandlboOll . iMiller fc Veilleuxl l2003h . Although it 
is well established that the power requirements to maintain these ionized layers can only 
be sup plied by stellar radiation, and not, for example, the mechanical input of supernovae 
shocks faevnoldslEggol ). this does not preclude the possibility of shocks altering the line ratios 
in the lower density, high z gas where the power requirements are much reduced. Shocks, 
of course, are just one example of non-photoionization heating, where the kinetic energy 
of gas is converted to thermal energy and radiation. Other possibilities include turbulence 
(kinetic energy — > thermal energy + radiation), magnetic reconnection (magnetic — )■ thermal 
energy + radiation), and cosmic-ray heating (particle kinetic energy — > thermal energy + 
radiation). 

In a stand ard plane-parallel shock ni odel, the optical line emission originates in three 
distinct zones (jPopita &: Sutherland! Il996l ). First, there is the non- equilibrium cooling zone. 
As the gas behind a shock front cools, it passes through the different ionization stages, where 
the ionization level is set by the gas temperature and by photoionization from the radiation 
produced by the hotter gas closer to the shock front. However, frequently the cooling is 
faster than the recombination and the ionization level is no longer in equilibrium. Ions with 
the highest ionization potential are the most likely to be produced in this zone. The second 
relevant region is the photoabsorption zone. For a strong shock, the gas is compressed by a 
factor of four by the shock front. As it cools, it compresses further, where the amount of 
compression depends on the relative importance of thermal and magnetic pressure. When the 
gas reaches ~ lO^K, it reaches a plateau in temperature until all of the ionizing flux produced 
in the post-shock gas has been absorbed. Since this zone is compressed, and irradiated by 
an extreme ultraviolet or X-ray flux, it usually dominates the optical line emission from 
the shock front, with emission line ratios distinctly different from standard H II regions. 
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For examp le. [Nelll] is enhanced relative to [QUI] by a factor of ~ 6 in shocks f Popita & 



Sutherland Il996l ) . The third region is the ionized precursor. Radiation from the post-shock 



gas streams in front of the shock, ionizing the low density pre-shock medium. 

There are two reasons why it is difficult to compare standard shock models to data on 
low density halos. First, these models are usually calculated in the absence of an external 
radiation field. Since the bulk of the emission for most of the optical/infrared lines tends 
to come from the photoabsorption zone (due to the higher emission measures), the halo 
radiation field must also be added self-consistently to the radiation emitted by the post- 
shock gas. For example, in the low density conditions of the halo, this photoabsorption zone 
may never fully recombine. 

The second difficulty is that at the low densities of the halo, the physical thickness of the 
shock increases to the point that the plane-parallel approximation becomes inappropriate. 
In can be shown that the emergent hydrogen ionizing photon fiux is proportional to the 
particle fiux through a plane-parallel shock front, = (p{vs)nH,psVs, where Vs is the shock 
speed, and UH^pg is the pre-shock density. The dimensionless efficiency factor, 4>{vs), depends 
principally on the shock speed and modestly on the amount of compression allowed by the 
magnetic field in the post-sh ock gas. This f actor ranges from 3 to 12 for shock speeds 



between 150 to 300 km s ^ (IBenjaminl 1 19941 ). This allows us to estimate the length of 



the precursor zone and the photoabsorption zone by assuming the ionizing photon fiux 
incident on each zone is balanced by the recombination in a gas column, I = $///(?T,|^a) = 
{(j)Vs) / {nHo) = (1200pc) {4> n-2)~^ VlooT^^'^ , where T4 is the gas temperature in units of 10^ 
K, fioo is the shock velocity in units of 100 km s~^, a = (2.6 x 10~^^ cm^ s~^)T4~°'^ is the 
hydrogen recombination rate, and n_2 is the hydrogen density in the irradiated gas in units 
of 10~^ cm~^. We see that for reasonable halo parameters, the thickness of the precursor 
zone can exceed the thickness of the entire ionized layer. For this reason, many shock model 
predictions exclude the contribution of the precursor. The thickness of the photoabsorption 
zone will be reduced by the compression factor experienced by the post-shock gas, typically 
by a factor of 4 to 400. The compression depends principally on the magnetic field in the 
post- shock gas. 

Because of the certain presence of the external radiation field, and the fact that the 
bulk of the optical line emission in shocks is reprocessed fiux from the post-shock gas, our 
previous approach of parametrizing the emission as a function of temperature, ionization 
parameter, and spectral hardness is still appropriate. However, the spectral shape will be 
much different than the blackbody flux that we have used. In principle, we could repeat the 
analysis in the previous subsection, replacing the blackbody spectrum with the spectrum 
emergent from shock-heated gas combined with the background stellar flux, and we will 
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consider this in future work. This would require the incorporation of (at least) two additional 
free parameters, the shock velocity/post-shock temperature and the flux level of ionizing 
photons. 



-39- 



4.3. Illustrative Physical Models for NGC 891 

In RWB we used the Monte Carlo photo-ionization code described in WM to attempt to 
understand whether the rise of the neon ratio with z in NGC 891 found in our Spitzer Cycle 
2 data could be explained by photo-ionization due to massive stars in a thin disk. We also 
considered the runs of the optical hne ratios shown in Figure 11. The details of the models are 
presented in RWB, therefore we emphasize here only a few of the most important features. 
Stars of temperature 40 kK, 45 kK and 50 kK were considered, with a total luminosity of 
10^^ ionizing photons per second, although one model considered a luminosity four times 
higher. While our standard models place the stars in a diffuse gas density structure (based 
on observations of NGC 891), one model also considers the effects of spectral hardening by 
denser gas around the stars - this involves passing the radiation field through a spherical gas 
distribution whose size is set to allow 15% of the ionizing photons to escape. 

Our main result was that no model could reproduce the rise of [Ne III] /[Nell] with 
z, with values of this ratio being generally overpredicted. The models did a better job in 
reproducing the observed runs of [SII]/IIa and [NII]/II«, but could not produce a rising 
[O III]/H/3 as observed, predicted too low values of [O I]/Hq;, and could not match the values 
of He I/Ha. 

The additional data on NGC 891 presented in this paper and the above exploration of 
parameter space in photo-ionization scenarios with CLOUDY show that the situation will 
not change - a thin disk of hot stars will not suffice. Instead they suggest that we vary the 
distribution of ionizing sources, such that gas at higher z is ionized primarily by sources that 
have, in the mean, harder spectra. Here we focus on such a source which is well suited to 
modeling by the Monte Carlo code by adding a secondary, hot stellar ionizing source with a 
vertically extended distribution. 

Our goal here is simply to present some illustrative models showing parameters for 
such a component (temperature, scale height, fraction of ionizing photons contributed) that, 
when added to our previous models, can approximately match the neon ratio behavior, 
rather than to perform an exhaustive search of parameter space. The first reason for this 
limited approach is that there are many possible ways in which the parameters suggested 
by our CLOUDY models could be incorporated into a physical model, not all of which 
may have a well-motivated physical explanation. Second, while we can find motivation for 
including the second source in several candidate types of stars based on Milky Way studies, 
information is scarce concerning their potential effectiveness as contributors to the ionization, 
and virtually non-existent regarding their distribution in NGC 891. Hence, our problem is 
not well constrained. 
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Hot sources in the Milky Way that may have a larger vertical extent than the main 
massive star layer include sdO and sdB stars, post-AGB stars, young white dwarfs, and 
runaway O stars. We return to the viability of each of these candidates after we present the 
modeling results. However, we do no t e here two related contexts in wh i ch su ch sources have 
been considered: iBinette et al.l (ll994r).lTaniguchi. Shiova. fc Muravamal (120001). and Stasihska 



et al. (120081 ) have explored possible importance of hot post-AGB stars and planetary nebula 
nuclei in explai n ing hi gh [0 III]/H/3 and jNH] /Ha ratios in Seyfert and LINER spectra, while 
Annibali et al.l ( l2010l ) and ISarzi et al.l ( 120101 ) have considered the importance of post-AGB 
stars for ionization in elliptical galaxies. 

Instead of the point-source geometry for the ionizing sources considered in RWB, here 
we use plane-parallel distributions. The massive stars occupy an infinitely thin layer, while 
an exponential describes the vertical distribution of the secondary component. We simulate 
a disk area of 4 kpc by 4 kpc. For the massive stars, we only consider a temperature of 30 
kK, and do not consider leakage. For the second component, we focus on hot, evolved, low 
mass stars and use stellar spectra from the Tiibingen NLTE model atmospheres fe.g. Ranch 



et al. I2OO9I ). We consider 50 kK models with elements from H to Ni, as well as ones with 
no H (He+C-|-N-|-0) to represent He-rich stars such as the sdOs. This choice should give 
us at least some idea of whether our result s are sensitiv e to composition, which can vary 
greatly among subdwarf stars, for example ( lHeberll2009l ). The total ionizing luminosity in 
our models is 5.5 x 10^° photons s~^. 

Figures 16 and 17 show resulting line ratios for two models where stars with the H- 
Ni spectra with a scale height of 1 kpc contribute 5% and 7.5% ("HNi5" and "HNi7.5" 
models) of the total ionizing luminosity, respectively. The extended component has only 
been included on the positive- 2; side, although it will also shghtly ionize the opposite side as 
the emission is isotropic. First, we note that almost all line ratios show a sharp spike or dip 
in the midplane which is not seen in the data. This is a consequence of the idealized nature 
of the models, where the midplane line of sight passes through the center of the ionized 
structure. Real lines of sight in the midplane presumably sample a range of such impact 
parameters, and should therefore reflect some mix of DIG and HII regions. Hence, we do 
not expect to match ratios correctly within the star forming layer. 



-41 - 




Fig. 16. — Vertical cuts showing the variation of hne ratios with z for the "HNi5" model (for 
doublets, the same line is plotted as in Figure 11). The vertically extended ionizing source is 
located on the positive z side only. The different curves in each panel represent cuts through 
the model centered at various in-plane coordinates, and reflect the noise in the simulation. 
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At z — 1 and 2 kpc, these models bracket the observed neon hne ratios seen in NGC 
891. Varying the contribution from the extended steUar component has httle effect on the 
behavior of [Nc II]/Ha with height, but greatly affects the behavior of [Ne II1]/Hq; (we do not 
report these ratios for the data mainly because extinction in the infrared and optical lines 
should vastly differ). Models with scale heights of 0.5 kpc and 2 kpc (not shown) produce, 
respectively, no rise, and much too steep a rise in the neon ratio off the plane. From this we 
can conclude that, for this representative temperature of the extended component, the rate 
of rise of the neon ratio between z — 1 and 2 kpc is governed mainly by the scale height of 
these stars, rather than their relative contribution to the ionizing flux. Hence, the only way 
to reproduce a rising neon ratio with such models is to include a vertically extended ionizing 
component, as we have done. 

Optical line ratios are more difficult to interpret, as already stated - nevertheless, we 
summarize here how these models fare in reproducing these ratios. Comparing Figures 11, 
16, and 17, the HNi5 model produces the rising trend in [0III]/H5 in the halo with values 
close to those observed, and comes reasonably close to matching the observed range of values 
of He 1/Ha of about 0.02 — 0.045. In these aspects, the model is much more successful than 
the ones presented in RWB. However, [S 11]/Hq; and [N11]/Hq;, although rising in the model, 
are now too low compared to the data. The models of RWB were more successful in this 
regard. [01] /Ha is also badly underpredicted, although this was also true for the RWB 
models. A negative abundance gradient with height will not help the situation: for instance, 
from the results of WM, [SIIJ/Ho; and [NIIJ/Ho; would be further lowered. 

The HNi7.5 model, while also producing rising [Sll]/Ha, [N11]/Hq; and [0111]/H/3 ra- 
tios, underpredicts the first two even more than the HNi5 model, and overpredicts the third. 
The behavior of He I/Ha is similar to that for the HNi5 model. 

The He-rich star model (not shown) produces a very similar behavior of the neon ratio 
to the HNi5 model when such stars contribute 30% of the ionizing flux. Their contribution 
is necessarily relatively large because the high helium content of the stars' atmospheres 
suppresses the neon- ionizing flux signiflcantly. However, this model yields very low [N II] /Ha 
(0.1 - 0.4) and [SII]/Ha (< 0.1) ratios that decline with z. [OIII]/H^ rises with z but is 
overpredicted, [O I] /Ha is badly underpredicted, while the behavior of He I/Ha is similar to 
that for the HNi5 and HNi7.5 models. We therefore do not consider such models further. 

If we focus only on the success of the HNi5 and HNi7.5 models in bracketing the observed 
behavior of the neon ratio with height, and ignore the problems introduced in matching other 
line ratios, we can then ask which type of aforementioned hot star could be responsible. The 
distribution of such stars is unknown in NGC 891, and therefore we can only try to constrain 
the possibilities based on what is known in the Milky Way. The only information available 
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on the distribution of stars in NGC 891 relative to the Milky Way is very general - thin and 
thick disk scale heights have been measured for both galaxies - and we end by summarizing 
these measurements. 

Of the various types of hot stars mentioned above, sdO stars may be the most likely 
in terms of temperature and z-distribution, but have the disadvantage that they tend to be 



He-rich. These stars have T = 40-80 kK, a. 



found to have T=45-55 kK (ITheill et al. 



;houg h in a sample of 21 He-rich sdOs, most were 



19941'). wh ile in another sample of 33 such stars, 



most were in the range 38-48 kK (jStroeer et al.ll2007l ). Temper atures for a sample of 13 sdOs 
with subsolar helium abundances are in the range 35 -80 kK (jStroeer et al.l 120071 ). Heights 
above the plane for the sample of iThejU et al.l (Il994j ) are estimated to be 0.4-1.5 kpc, and 
ThejU et al . specu late that they may form a thick disk or halo population. More recently. 



Napiwotzkil ( 120081 ) finds from kinematics of 31 He-sdOs that 41% are in the thin disk, 39% 



in the thick disk, and 20% in the halo. Their space density is very uncertain. iDownesI (119861 ) 
found a value in the Solar neighborhood of about 7 x 10^^ pc~^, although with a low scale 
height of 100 — 200 pc. If the scale height is indeed 1 kpc, the midplane space density would 
be reduced. The midplane space density of such stars in our HNi5 model can be derived 
if we assume that our extend ed componeri t cons ists of stars of a typical sdO luminosity of 



150 L0, based on values from lThejU et al.l ( 119941 ). We find a value of 6 x 10 ^ pc ^. Hence 



it seems, although information is very scarce, that such an extended component cannot be 

ruled out on grounds of space density, temperature or scale height. 
— I 



de Boer et al. ( 119971 ) find that the sdB stars also have a large scale height of about 1 



kpc, but these stars are certainly too cool to be responsible f o r the neon line behavior. The 



Milky Way post-AGB stars in the catalog of ISzczerba et al.l (120071 ) are also too cool, with 



(Green. Schmidt, & '. 


jiebert 


1986 


), which have temperatures of order 100 kK. Only 40 


known ( 


Werner et al. 


2008 


)• 



A sample of white dwarfs in the Sloan Digital Sky Survey are fo und to have a scal e 
height of about 340 pc, consistent with several earlier determinations (iHarris et al.ll2006l ). 
These authors also find a possible trend of increasing scale height with decreasing luminosity 
(and therefore temperature), but the statistical significance is uncertain. Such a trend would 
naturally be expected because of heating of the disk over time. Nevertheless, it seems that 
luminous, hot white dwarfs should have a scale height of much less than 1 kpc. If the PG1159 
stars are evolving into hot white dwarfs, the same conclusion should apply to them. 



Runaway O stars account for 10 — 30% of all O stars ( lGieslll987l ) and form a distribution 
more vertically extended than the thin disk. However, most have V/s^'s of less than 100 km 
s~^, so it is difficult to see how they could form a layer with a 1 kpc scale height given their 
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lifetimes. In fact, such stars in the sample of iMoffat et al.l (119981) are mostly within 200 
pc of the midplane. Some B stars are found at greater heights (lAUen fc KinmanI l2004j . and 
references therein), but again are too cool to be candidates here. 

Of course, candidate sources with appropriate temperatures may have a different scale 
height and volume density in NGC 891. While individual candidate stellar types are difficult 
to study, the o verall stellar disk struct u re of NGC 891 has been modeled several times, 
beginning with Ivan der Kruit fc Searld (119811). w ho were the first to report a thick disk. 



From i?-band photometry, iMorrison et al.l (jl997l ) derived scale heights for the thin and 
thick disks of 400 — 650 pc and 1.5 — 2.5 kpc, resp ectively, with the thi ck disk contributing 
5 — 10% of the surface brightness in the midplane. IXilouris et al.l ( 119981 ). modeling starlight 
and exti nction from five-band photona e try, f ound a thin disk scale height of 400 pc. More 
recently, llbata. Mouhcine. &: Rejkubal ( 120091 ) found a thick disk scale height of 1.44 ± 0.03 
kpc from star counts in HST ACS ir nages. Thin and t hick disk scale heights are therefore 
larger than in the Milky Way, where lJuric et al.l ( l2008l ) find values of 0.3 kpc and 0.9 kpc 
from the Sloan Digital Sky Survey. Hence, there may well be an extended distribution of 
hot, evolved stars in NGC 891, but one cannot currently say more than this. 

We conclude from this initial exploration that such a component cannot be ruled out as 
an explanation for the neon ratio behavior given what is known about the stellar distribution 
in NGC 891 and the required space density of this component in our models. It also may 
help explain the [O III]/H/3 and He I/Ha behavior in this galaxy. However its inclusion does 
introduce problems with other optical line ratios. As previously mentioned, a more complete 
exploration of parameter space in such models, including a consideration of non-ionizing 
heating, may result in a better match to the data for NGC 891, but is beyond the scope of 
this paper. 

Finally, we note that any contribution of a hidden AGN or nuclear super star cluster 
with a hard spectrum is very unlikely to be a factor in NGC 891. First, the structure of the 
Ha halo emission correlates with star formation in the disk, both in the contrast on either 



side of the minor axis, and on the scale of individual filaments (lRandlll998al ). Second, in a 
slit running parallel to the disk at a height of 700 pc in the DIG halo, [NII]/Ha, [S II]/Ha, 
and [01] /Ha all anti-correlate with Ha halo intensity, dropping at the location of filaments 
(jRandlll998bl ). These results imply a strong connection between line ratios and disk star 
formation, with no indication of a contribution from a nuclear source. 



- 46 - 



Discussion of PAHs 



Our main results for PAH emission scale heights are as follows. First, the scale heights 
(depending on extinction) are many hundreds of pc. This res ult adds to the growing e vidence 
for v ertically extended PAH distributions in galaxies fe.g.. lEngelbracht et al.ll2006l: Beirao 
et al. 



2008 



Howt 



d l2009h ■ 



Spec ifically for these two galaxies, in the IRAC 8 fj,m image of 



NGC 891, IWhaley et al.l ( l2009l ) measure a scale height of 248 ± 50 pc (with no extinction 
correction), while the IRAC 8 /xm image of NGC 5775 presented by lHowkl (120091 ) clearly shows 
a very extended layer with much vertical filamentary structure - although this wavelength 
traces more ionized PAHs. Second, the PAH scale heights are greater in NGC 5775 than in 
NGC 891 - the former also fe aturing more active star formation and more extended warm 

Rand 



ioniz ed and neutral gas halos (llrwin 



1994 



Collins et al. 



2000 



Oosterloo et al. 



2007 



et al. Il990l ). Our PAH scale heights are comparable to or somewhat lower than the HI scale 
heights. Third, within NGC 891 and NGC 5775, scale heights are similar among the various 
PAH features. 

There are two obvious caveats to consider when interpreting these PAH scale heights. 
The first is that we have halo measurements at only one distance along the major axis in 
each galaxy. The second is that the scale heights should refiect the intensity of the exciting 
radiation as well as the PAH density distribution. 

No PAH feature stands out as having a significantly greater or smaller scale height in 
either galaxy halo. Thus, there is little evidence for significant modification between disk 
and halo of the relatively large, neutral PAH population responsible for the discrete features 
in this wavelength range. However, for NGC 891, the 8 Atm (dominated by the PAH feature 
at 7.7 fim) scale height of 248 ± 50 pc measured by IWhaley et al.l ( 120091 ) (uncorrected for 
extinction) is significantly lower than our average of 475 ± 30 pc for the PAH features. In 
order for this difference to be an extinction effect, the extinction at 8 /im would need to be 
about 1 mag greater than for our PAH features. Although 9.7 /xm silicate absorption may 
have a significant effect in the 8 /zm band, it may also affect the 11.2 p, m band (perhaps eve n 



more so than the 7.7 fim feature which dominates the 8 fim images; iGalliano et al. 



l2008h . 



while the 17.4 /xm feature could suffer from significant, albeit less, extinction due to the 18 
Hm silicate feature. We therefore consider extinction an unlikely explanation. Furthermore, 
the PAH emission spectrum should be relatively independent of ISRF strength for the weak 



radiation fields ex pected in integrated disk and halo lines of sight (IGalliano et al.l 12008 



Draine fc Lill2007l ). Therefore the difference in scale heights most likely represents a real 
change in the PAH population between disk and halo. The result seer ns consistent with th e 
drop of the 7.7 /im/11.2 fim intensity rati o with z in the M82 outflow (IGalliano et al.ll2008l ). 
According to the models of iDraine fc Lil (120071 ) . the 7.7 /im feature is due to smaller and 
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more ionized PAHs than those responsible for the 10 — 20 um features, although Galliano 



et al. ( 120081 ) claim that the 7.7 /im/11.2 fim ratio is sensitive only to ionization and not 
size variations. Our results would then indicate a significant drop in PAH ionization, and 
possibly a transition to larger PAHs, with height in the halo of NGC 891. 

These results suggest a picture in which PAHs are transported upward from the disk 
(where they are thought to form in C-rich evolved AGB star envelopes) coupled with warm 
gas as part of a disk-halo fiow. However, as disc ussed in ^1, o t her rn echanisms could be at 
work. Regarding radiative acceleration of dust, iDavies et al.l ( 1l998l ) find that the process 
is less efficient for smaller grains. Alternatively, PAHs could form in situ. One possible 
mechanism is the shatterin g of larger grains in shocks, which could be effective for shock 



speeds as low as 50 km s ^ ( jjones et al 



119961 



i 



Shocks may also destroy PAH s, with larger ones being able to withstand faster shocks. 
Micelotta. Jones, fc Tielend (120101 ) find partial destruction of PAHs for shocks of 75 — 100 
km s~^, while ones with 50 and 200 C atoms do not survive shocks faster than 100 km s~^ 
and 125 km s~^, respectively. PAHS in the 10 — 20 fim range probably contain a few hundred 
to a few thousand C atoms (jPraine fc Lil 120071 ). so if shock speeds are less than about 125 
km s~^ they should survive. On the other hand, the lower scale height of the 8 /im emission 
relative to the 10 — 20 fim features co uld be due to shock modification. If the former arise 
mostly from PAHs of ~ 100 C atoms ( iDraine &: Lil 120071 ) there ma y be partial destruction 
for shock speeds near 100 km s"^, if one interpolates the models of iMicelotta et al.l (120101 ) 
for the 50 and 200 C atom cases. 

We have also found that most 10 — 20 PAH EWs are higher at z = 1 — 2 kpc t han in 



the midplane, using two different methods to estimate the continuum. iBeirao et al.l (120081 ) 
find a similar trend between disk and halo for PAH EWs in this wavelength range in M82. 
One might conclude that, together with the increase in [Ne III]/[Ne II] seen in NGC 891 and 
NGC 5775, th e result reflects p revious flndings of an increase of EWs with radiatiori fleld 
hardness (e.g., IWu et al.ll2006l for Blue Compact Dwarf Galaxies; lEngelbracht et al.l 12008 



le-g. 

for starbursts 

correlation does not exist in star forming regions (IBrandl et al.ll2006 



However, our values are in the r egime [NeIII]/[Ne 



PAH features also weaken in galaxies with low metallicitv (e.g 
et al. 



I j < 1, where such a 



Gordon et al.l 120081) . 



Madden et al.l 



20061: Wu 



20061 ) . There is no direct information on how metallicity varies with height in our 



three halos. Some line s of evidence, for instance the generally low metallicities of Milky Way 
High Velocity Clouds (IWakkerl 120041 ). suggest halos may have lower metallicities than disks, 
possibly as a result of mixing of disk-halo cycled a nd primordial gas. H owever, these clouds 
are found at heights much greater than 1 — 2 kpc (IWakker et al.ll2007l ). Regardless, a drop 
in metallicity with height would be expected to produce, if anything, a drop in EWs. 
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Hence, the explanation of t his result is not clea r. As mentioned before, it should be kept 
in mind that there is evidence ( iDraine et al.l 120071 ) from SINGS galaxies that, just like the 
discrete features, the continuum in the 16 /im range is mainly due to single-photon heating 
of PAHs. In fact, feature and continuum radiation in this range is dominated by molec ules 
of size 5 — 20 A molecules for weak radiation fields in the model of iDraine fc Lil (120071 ) . If 
the continuum throughout the 10 — 20 /xm range is also dominated by single-photon PAH 
heating, then our result suggests a perhaps rather subtle change in the PAH population 
between disk and halo. It is not clear which, if any, of the physical mechanisms described 
above may be responsible. 



6. Conclusions 

We have presented 10 — 20 /im spectroscopy of the disks and halos of the edge-on galaxies 
NGC 891, NGC 5775 and NGC 3044. Regarding the gas-phase lines, our main result is that 
[Ne III] /[Nell] is higher in the halos of NGC 891 and NGC 5775 than in the disks. Scatter 
in NGC 3044 prevents any trend from being seen. 

These results exacerbate the problem of explaining DIG line ratio behavior with simple 
photo-ionization models featuring a thin disk of massive stars. Focusing on NGC 891, where 
we have the most observational constraints, we have explored parameter space in CLOUDY 
photo-ionization models to try to determine, as a function of height, what combination 
of radiation temperature (assuming blackbodies for simplicity), ionization parameter and 
gas temperature may reproduce the observed IR and optical line ratios. We find that a 
dramatic rise in radiation temperature with z is required, as well as a more modest rise in 
gas temperature and a fall in ionization parameter. Not all line ratios can be reproduced 
without considering more complex models, however. We then considered representative 
physical photo-ionization models incorporating a thin disk of massive stars in combination 
with a thick disk of hot (50 kK), presumably evolved stars. We found that the neon ratio 
behavior in particular could be approximately reproduced for such a thick disk with scale 
height of order 1 kpc and a contribution of 5 — 7.5% of the ionizing radiation. Whether 
such a component exists in NGC 891 (or even the Milky Way) remains very unclear, but 
it cannot yet be ruled out based on space density, temperature and scale height arguments. 
In our models, such a component is also reasonably successful in reproducing the rising 
[0III]/H/3 and approximate He I/Ha values in this galaxy. However its inclusion does 
introduce problems with other line ratios. A secondary stellar source is not the only way 
of increasing radiation hardness with height, and our exploration of parameter space with 
CLOUDY suggests a feasible way in which other such sources, including hard EUV/X-ray 
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radiation, may be included in future work. 

PAHs in NGC 891 and NGC 5775 form a vertically extended distribution, with emission 
scale heights larger in NGC 5775, as is the case for the HI gas. The PAH emission scale 
heights are comparable to, or somewhat less than, the HI values. These results suggest 
that PAHs participate in the active disk-halo flows in these two galaxies. In NGC 891, the 
scale heights are larger than the 8 /xm scale height, suggesting a transition to more neutral, 
and possibly larger, grains with z. Most PAH equivalent widths are higher in the halos 
than in the disks, again indicating some slight modification in the PAH population in the 
halo environment. Shocks and radiative acceleration are processes that may affect PAH 
populations in halos. 

The 17 [xm H2 line is detected in all pointings, indicating a warm molecular gas compo- 
nent with a large vertical extent, reaching as high as 2; = 2 kpc. The greater scale height of 
this emission in NGC 5775 vs. NGC 891 again suggests a connection with disk-halo flows. 
Column densities and temperatures cannot be derived without measurement of additional 
transitions, however. Future observations will no doubt reveal the fraction of halo gas mass 
in this phase, and its importance for understanding the ISM of halos. 

Finally, we wish to stress here the importance of the neon ratio in understanding the 
energetics of the diffuse ISM, and the great advantage it holds over optical diagnostics. 
Future instruments capable of sensitive, high spatial and spectral resolution observations of 
these and other mid-IR lines should be a strong priority as these fines reveal a wealth of 
information. 
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Table 1. Emission Lines and Ionization Potentials 



Line 


Ionization Potential (eV) 


Ha 6563A 


13.6 


[SII] 6716A 


10.4 


[Nil] 6583A 


14.5 


[OI] 6300A 


13.6 


[OIII] 5007A 


35.1 


He I 5876A 


24.6 


[NeH] 12.81 fj,m 


21.6 


[NeHI] 15.56 //m 


41.0 


[SHI] 18.71 nm. 


19.1 


[SIV] 10.51 nm 


34.8 
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Table 2. IRS Observations 



Date 


Pointing'' 


R.A. (J2000.0) 


Decl. (J2000.0) 


Integ. time per nod (sec) 


2008 Mar 26 


NGC 5775 Disk 1 


li^ 53™ 58.4= 


3° 32' 24" 


975 


2008 Mar 26 


NGC 5775 Disk 2 


14^ 53™ 59.0= 


3° 32' 12" 


975 


2008 Mar 26 


NGC 5775 Disk 3 


14!^ 53™ 59.6= 


3° 32' 00" 


975 


2008 Mar 26 


NGC 5775 Halo 


14'' 53™ 57.5= 


30 i^// 


3854 


2008 Mar 26 


NGC 5775 Sky 


14^ 54™ 11.7^ 


3° 35' 58" 


3854 


2008 Sep 9 


NGC 891 Halo east 


2^ 22™ 40.5= 


42° 22' 13" 


21194 


2008 Sep 9 


NGC 891 Halo west 


2^ 22™ 33.2= 


42° 22' 44" 


21194 


2008 Sep 9 


NGC 891 Sky 


2^ 23™ 32.7= 


42° 22' 21" 


21194 


2009 Jan 7 


NGC 3044 Halo BE 


9^ 53™ 39.1= 


1° 34' 42" 


10597 


2009 Jan 7 


NGC 3044 Disk 1 


9^ 53™ 38.6= 


1° 35' 01" 


975 


2009 Jan 8 


NGC 3044 Disk 2 


9'' 53™ 39.5= 


1° 34' 55" 


975 


2009 Jan 8 


NGC 3044 Disk 3 


9'' 53™ 40.4= 


1° 34' 49" 


975 


2009 Jan 8 


NGC 3044 Halo NW 


9'' 53™ 39.9= 


1° 35' 08" 


10597 


2009 Jan 8 


NGC 3044 Sky 


9'' 53™ 39.5= 


1° 36' 30" 


10597 



Pointing center of the field of view (each nod pointing is symmetrically offset from this position) 
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Table 3. NGC 891 Infrared Line Intensities and PAH Equivalent Widths 





East Halo 


East Halo 


West Halo 


West Halo 


Line 


Intensity^ 


EW 


Intensity 


EW 












[SIV] 10.51 nm 


<0.17'' 




<0.16 




[Nell] 12.81 fim 


2.4±0.2 




3.0±0.2 




[NcIII] 15.56 Atm 


1.6±0.1 




1.4±0.1 




[SIII] 18.71 Atm 


1.7±0.1 




l.OitO.l 




H2 5(1) J = 3 - 1 17.03 Atm 


1.4±0.1 




1.7±0.1 




PAH 11.2 nm 


18.2±1.3 


1.14±0.04 


19.2±1.4 


1.36±0.05 


PAH 12.0 nm 


1.6±0.3 


0.10±0.02 


2.2±0.3 


0.15±0.02 


PAH 12.7 nm 


7.1±0.5 


0.52±0.02 


7.8±0.6 


0.64±0.02 


PAH 16.5 fan 


O.SitO.l 


0.11±0.02 


1.2±0.1 


0.23±0.02 


PAH 17.4 nm 


0.4±0.1 


0.07±0.02 


0.7±0.1 


0.13±0.02 


PAH 10.8 - 14.5 Atm 




4.45±0.13 




0.68±0.02 


PAH 16.3 - 17.7 fim 




4.45±0.13 




0.68±0.02 


15 /im continuum'^ 


7.0±0.5 




0.6.3±0.4 





''Intensity units are lO"^'^ erg cm~^ arcsec"^ 
''Upper limits are 3a 
''Integrated over 14.5 — 15.5 /nm 



Tabic 4. 


NGC 891 Line Ratios 


Position 


[Ne III] /[Nell] 


[SIII] /[Nell] 


East Halo 


0.67±0.04^ 


0.71±0.04 


West Halo 


0.47±0.03 


0.33±0.03 


'^Errors are 


based on random 


noise in the data 



and the wavelength-dependent pointing uncer- 
tainty only. 



Table 5. NGC 5775 Infrared Line Intensities and PAH Equivalent Widths 





Disk 1 


Disk 1 


Disk 2 


Disk 2 


Disk 3 


Disk 3 


Halo 


Halo 


Line 


Intensity" 


EW 


Intensity 


EW 


Intensity 


EW 


Intensity 


EW 






(/im) 




(Aim) 




(Aim) 




(Aim) 


[SIV] 10.51 fim 


<3.7'' 




<3.8 




<3.9 




<0.5 




[NcII] 12.81 ^tm 


96±7 




127±8 




155±10 




13.8±0.9 




[Nelll] 15.56 /xm 


17±2 




24it2 




28±2 




3.2±0.3 




[Bill] 18.71 nm 


31±2 




39±3 




56±4 




4.3±0.3 




H2 5(1) J = 3 - 1 17.03 ixm 


29±2 




32±2 




30±2 




3.0±0.2 




PAH 11.2 fim 


1045±72 


1.05±0.03 


1006±70 


1.03±0.03 


972±67 


1.05±0.03 


84±6 


1.12±0.04 


PAH 12.0 fim 


60±7 


0.05±0.005 


64±7 


0.05±0.005 


70±8 


0.06±0.006 


10±1 


0.12±0.009 


PAH 12.7 Atm 


356±24 


0.42±0.01 


338±23 


0.40±0.01 


368±25 


0.51±0.01 


23±2 


0.41±0.01 


PAH 16.5 Aim 


55±4 


0.11±0.007 


57±5 


0.12±0.008 


54±5 


0.12±0.008 


4.4±0.4 


0.11±0.009 


PAH 17.4 urn. 


28±3 


0.06±0.005 


23±3 


0.04±0.005 


21±3 


0.04±0.006 


3.9±0.4 


0.10±0.009 


PAH 10.8 - 14.5 Aim 




3.4±0.1 




3.6±0.1 




3.75±0.1 




3.83±0.1 


PAH 16.3 - 17.7 ixm 




0.68±0.02 




0.78±0.02 




0.58±0.01 




0.71±0.02 


15 Aim continuum'^ 


427±28 




427±28 




402±26 




33±2 





"Units arc lO"^'^ erg cm~^ arcsec"^ 
''Upper limits arc 3a 
'^Integrated over 14.5 — 15.5 A^m 
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Table 6. NGC 5775 Line Ratios 



Position 


[Ne III] /[Nell] 


[S III] /[Nell] 


Disk 1 


0.18±0.02=' 


0.32±0.01 


Disk 2 


0.19±0.01 


0.31±0.01 


Disk 3 


0.18±0.01 


0.36±0.01 


Halo 


0.23±0.01 


0.31±0.01 



"Errors are based on random noise in the data 
and the wavelength-dependent pointing uncer- 
tainty only. 



Table 7. NGC 3044 Disk Infrared Line Intensities and PAH Equivalent Widths 





Disk 1 


Disk 1 


Disk 2 


Disk 2 


Disk 3 


Disk 3 


Line 


Intensity" 


EW 


Intensity 


EW 


Intensity 


EW 






(/im) 




(/im) 




(/im) 


[SIV] 10.51 fim 


<1.1 




<2.4'' 




<4.3 




[Nell] 12.81 Atm 


30.8±2.1 




82±5 




184±12 




[Nelll] 15.56 urn 


18.5±1.2 




61±4 




54±4 




[Bill] 18.71 nia 


11.9±0.8 




27±2 




48±3 




H2 5(1) J = 3 - 1 17.03 /im 


11.2±0.8 




12±1 




38±3 




PAH 11.2 (um 


201±14 


1.05±0.03 


561±39 


0.82±0.02 


1066±74 


0.84±0.02 


PAH 12.0 Atm 


31±3 


0.17±0.01 


33±5 


0.05±0.01 


77±8 


0.05±0.01 


PAH 12.7 Atm 


82±6 


0.68±0.01 


215±15 


0.46±0.01 


423±27 


0.44dz0.01 


PAH 16.5 fj.m. 


7.2±0.9 


0.09d=0.01 


12±2 


0.04±0.01 


73±6 


0.12±0.01 


PAH 17.4 /im 


4.6±0.8 


0.06±0.01 


15±2 


0.06±0.01 


23±3 


0.04±0.01 


PAH 10.8 - 14.5 nm 




3.2±0.1 




2.8±0.1 




3.1±0.1 


PAH 16.3 - 17.7 ixm 




0.42±0.01 




0.23±0.07 




0.49±0.01 


15 /um continuum'^ 


75.3±4.9 




276±18 




527±34 





''Units are lO"^'^ erg cm"^ s~^ arcsec"^ 
''Upper limits arc 3cr 
"^Integrated over 14.5 — 15.5 /im 
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Table 8. NGC 3044 Halo Infrared Line Intensities and PAH Equivalent Widths 





North Halo 


North Halo 


South Halo 


South Halo 


Line 


Intensity"' 


EW 


Intensity 


EW 






/im 




/xm 


[SIV] 10.51 urn 


<0.18 




<0.18'' 




[NcII] 12.81 fim 


4.8±0.3 




4.4±0.3 




[Nc III] 15.56 fim 


1.7±0.1 




3.7±0.2 




[SIII] 18.71 /im 


2.4±0.2 




1.23±0.09 




H2 S{1) J = 3 - 1 17.03 urn 


3.7±0.2 




1.9±0.1 




PAH 11.2 /im 


36.9±2.6 


1.15±0.04 


33.5±2.3 


1.28±0.04 


PAH 12.0 /im 


2.7±0.3 


0.08±0.008 


6.1±0.5 


0.22±0.01 


PAH 12.7 fj,m 


9.6±0.7 


0.40±0.01 


12.5±0.9 


0.57±0.02 


PAH 16.5 fim 


3.8±0.3 


0.21±0.01 


2.4±0.2 


0.15±0.01 


PAH 17.4 fim 


1.2±0.2 


0.07±0.01 


2.5±0.2 


0.17±0.01 


PAH 10.8 - 14.5 fim 




3.4±0.1 




3.3±0.1 


PAH 16.3 - 17.7 urn 




0.74±0.02 




0.68±0.02 


15 jUm continuum'^ 


IS.litl.O 




12.6±0.8 





^Units arc 10^"'^'' erg cm~^ arcsec"^ 
''Upper limits are 3cr 
■^Integrated over 14.5 — 15.5 /im 



Table 9. NGC 3044 Line Ratios 



Position 


[Ne III] /[Nell] 


[SIII] /[Nell] 


Disk 1 


0.60±0.02» 


0.39±0.01 


Disk 2 


0.74±0.02 


0.32±0.01 


Disk 3 


0.29±0.01 


0.26±0.01 


North Halo 


0.35±0.02 


0.49±0.02 


South Halo 


0.84±0.03 


0.28±0.01 



■''Errors are based on random noise in the data 
and the wavelength-dependent pointing uncer- 
tainty only. 
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Table 10. Equivalent Exponential Scale Heights for Gas and PAHs in NGC 891 and NGC 

5775, Uncorrected for Extinction 



Line 


NGC 891 East Halo 


NGC 891 West Halo 


NGC 5775 


[Nell] 12.81 At 


420 ± 100 


400 ± 50 


890 ± 180 


[Ne III] 15.56 n 


640 it 190 


540 ± 80 


1000 ± 260 


[SIII] 18.71 n 


530 ± 240 


470 ± 40 


860 ± 190 


H2 S{1) J = 3 - 1 17.03 nm 


580 ± 160 


550 ± 80 


850 ± 60 


PAH 11.2 ^lT[^ 


470 ± 60 


430 ± 10 


790 ± 50 


PAH 12.0 


530 ± 250 


510 ± 80 


1060 ± 150 


PAH 12.7 Atm 


480 ± 80 


440 ± 40 


720 ± 80 


PAH 16.5 fim 


460 ± 40 


460 ± 20 


780 ± 60 


PAH 17.4 nm 


470 ± 30 


500 ± 50 


1080 ± 220 


PAH 15 //m continuum 


440 ± 100 


380 ± 20 


780 ± 40 



Table 11. Derived physical parameters vs. height in the halo of NGC 891 



Height 


log!7 


logT 


logTi, 






T 


n 


(kpc) 








(10*^ ph s-l cm- 




(K) 


(K) 


z = 0'' 


-3.10 ±0.10 


4.56 ± 0.03 


3.85 ± 0.05 


(2.40 ± 0.55)nH,. 


-1 


7130 ± 820 


36, 000 ± 2000 


z = l 


-3.48 ± 0.06 


4.66 ± 0.03 


3.86 ± 0.03 


(1.00±0.14)nH^,. 


-1 


7260 ± 500 


45, 800 ± 3200 


z = 2 


-3.80 ±0.05 


4.76 ± 0.03 


3.94 ± 0.02 


(0.48 ± 0.05)njf,. 


-1 


8600 ± 300 


57, 000 ± 3300 



^As derived from matching [Neni]/[Nen], [SII]/Ha, [NHj/Ha, and [Oni]/Ha. The difficulties of matching [OI/Ha, Hel 
5876/Ha, and [Sni]/[Nen], and implications thereof, are discussed in the text. 

''All fluxes of hydrogen ionizing photons are scaled to the same density. Assuming a gas scaleheight of hz = 2.3 kpc and no 
change in gas flUing fa<;tor would reduce the fluxes at z = 1 kpc and z = 2 kpc by a factor of 0.65 and 0.42 respectively. 

'^Emission line ratios come closest together at listed parameters, but there is no unique solution. 



